arXiv: 1501.04364v2 [nlin.CD] 28 Feb 2015 


Thirty Years of Turnstiles and Transport 

J. D. Meisf] 

Department of Applied Mathematics 
University of Colorado 
Boulder, CO 80309-0526 

(Dated: March 3, 2015) 


To characterize transport in a deterministic dynamical system is to compute exit time distributions 
from regions or transition time distributions between regions in phase space. This paper surveys 
the considerable progress on this problem over the past thirty years. Primary measures of transport 
for volume-preserving maps include the exiting and incoming fluxes to a region. For area-preserving 
maps, transport is impeded by curves formed from invariant manifolds that form partial barriers, 
e.g., stable and unstable manifolds bounding a resonance zone or cantori, the remnants of destroyed 
invariant tori. When the map is exact volume preserving, a Lagrangian differential form can be 
used to reduce the computation of fluxes to finding a difference between the action of certain key 
orbits, such as homoclinic orbits to a saddle or to a cantorus. Given a partition of phase space into 
regions bounded by partial barriers, a Markov tree model of transport explains key observations, 
such as the algebraic decay of exit and recurrence distributions. 

PACS numbers: 05.45.-a, 45.20.Jj, 47.52.+j 


The problem of transport in dynamical systems is to quantify the motion of collections of trajec¬ 
tories between regions of phase space with physical significance, for example, to determine chemical 
reaction rates, mixing rates in a fluid, or particle confinement times in an accelerator or fusion 
plasma device. For Hamiltonian systems with a phase space containing regular (elliptic islands and 
tori) and irregular (chaotic) components, transport is impeded by partial barriers bounding reso¬ 
nance zones or by remnants of invariant tori. In the former case the barriers are formed from the 
broken separatrices of an unstable periodic orbit, and in the latter, they are formed from similar 
stable and unstable manifolds of a remnant torus—a cantorus. Thirty years ago, MacKay, Meiss, 
and Percival discovered that cantori form robust partial barriers for area-preserving maps: in any 
region of phase space the most resistant barriers are the cantori. In this review, I discuss this 
work and survey subsequent developments. While much has been done, there are still many open 
questions. 


I. INTRODUCTION 

In this article, I attempt to survey the history of the theory of transport in deterministic, conservative dynamics, 
as well as to discuss some of the many questions that remain. My interest in this subject arose from a collaboration 
with Robert MacKay and Ian Percival that began just over thirty years ago with the publication of two papers on 
“Transport in Hamiltonian Systems” |MMP84al IMMP84b1 . Like many fruitful ideas, there was some synchronicity 
to this discovery: we learned at the “Dynamics Days” meeting in June of 1983 1 that David Bensimon and Leo 
Kadanoff were working on a similar approach [BK84j . I remember that when Robert, Ian, and I celebrated the 
acceptance of our papers, we thought that they contained ideas for at least ten years of additional research. Some 
thirty years later, I believe there are still many fruitful, open problems in this held. Much has been done: our 


* Electronic address: jdm@colorado.edu 

1 This was the sixth Dynamics Days and the fourth in Enschede, organized Robert Helleman. 
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articles have been cited 697 times, according to the Web of Science. The work of Vered Rom-Kedar and Stephen 
Wiggins, which introduced the term lobe dynamics and (with Anthony Leonard) applied these ideas to fluid mixing, 
[RKW9011RKLW90] has also been extensively cited (443 times). 

The problem of deterministic transport also has been studied frequently in the first 24 volumes of Chaos: An 
Interdisciplinary Journal of Nonlinear Science. Indeed the terms turnstile or lobe have been mentioned in 114 
articles, beginning with an article on transport in a four-dimensional map in the very first issue [ASZ91I . It would 
be impractical to try to review all of this work. Since I wrote a review in 1992 |Mei92| . I will concentrate here 
on describing some of the more recent advances, as well as on posing some questions that I hope will be answered 
before the 50 th anniversary of Chaos. However, to make this paper coherent, it is necessary to tell some of the 
early story. 

A theory of transport seeks to characterize the motion of groups of trajectories from one region of phase space to 
another. Typical dynamical systems have both regular and chaotic (irregular) trajectories. The question “Where 
does this trajectory go?” is a particularly difficult one to answer for chaotic trajectories since they—by definition- 
exhibit sensitive dependence upon initial conditions. The intermixture of regular and irregular dynamics is par¬ 
ticularly evident in the Hamiltonian case where stable periodic orbits are surrounded by elliptic islands that are 
embedded in the chaotic sea as a fractal. 

So how does one study transport? Our original paper lMMP84bi set out a method: 

In this paper we study transport in the irregular components and initiate a theory of the organization 
inherent in the apparent chaos. This is achieved by recognizing a natural division of the irregular 
components into regions separated by partial barriers formed by joining the gaps in invariant Cantor 
sets. We discuss the mechanism of transport between these regions. 

There were many antecedents to these ideas. For example, in his paper introducing the quadratic area-preserving 
map, Henon computed the set of orbits that do not escape to infinity, a kind of anti-transport result |Hen69j . More 
direct antecedents include the development of transition state theory in chemistry [PP791 IWW10I IMS 14 and 
especially the notions of Wigner for finding surfaces of minimal flux |Wig371 IKec67j . For a chemical reaction that 
has no activation energy and for which classical mechanics is a good approximation, Wigner thought there should 
be a surface in phase space that distinguishes between dissociated and reacted states. Unless the surface has 
“wrinkles”, that is, unless there are trajectories that cross it twice during one collision (similar to what MacKay 
later called “sneaky returns” |Mac94 | ). then the one-way flux across this surface gives a bound on the reaction rate. 

For the purposes of this exposition consider a map / : M —► M on some n-dimensional manifold M. Most of the 
ideas here can be applied directly to flows, either directly using a Poincare or stroboscopic section, or by extension 
to continuous time |MM861 [Mac941 IMSM13I . For a map, a trajectory is a sequence {z t } = {..., Zj,Zj+ 1 ,...}, such 
that Zt £ M and 


z t+ i =f(z t ). 

We will always assume that / is a homeomorphism, and typically that it is a diffeomorphism. In addition, we will 
suppose that / preserves a volume measure /i: for any measurable set A, 

M/0 4 )) = MM>- (!) 

When the map is a diffeomorphism, volume can be computed using a differential form 

= p{z)dz x A dz 2 A ... A dz n , (2) 

with density p 1 i.e., p{A) = f A H. Then / preserves volume when 

f*n = n. (3) 

Here the pullback , /*, is the local action of the map / on differentials. An easy mnemonic for how this works 
is to denote the image as a function z'(z) = f(z), and then f*dz = dz' is the differential of this function, i.e., 
f*dz l = D .h.j{z)dzC Here Dfi.j = df l jdzj is the Jacobian matrix of /. The implication is that |3j) for the 
form ([2]) implies that 

det (Df{z))p{f{z)) = p(z). 


For the common case that p{z) = 1, this implies that the Jacobian of / has determinant one. 
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II. AREA-PRESERVING MAPS 

Most of the theory of transport has been developed for the area-preserving case. For an area-preserving map , 
there are coordinates z = (cc, z/) such that the preserved volume form is tt = dy A dx. 2 The well-studied Chirikov 
standard map on the cylinder M = Sxl, 

(x 1 , y') = (x + y' mod 1, y - £ sin(27ra)) , (4) 

is one such family. A typical phase portrait for this map is shown in Fig. |Tj The standard map is an example of 
a twist map on the cylinder: a map for which the image of a vertical line is a graph over the angle |Mei92| . The 
powerful Aubry-Mather theory applies to this class of maps jMat82l IAub83| . Another example is Henon’s family 
of area-preserving maps, 

{x’,y') = {~y + 2(a - x 2 ),y). (5) 



FIG. 1: Phase portrait of the standard map 0 for k = 0.975, on the domain [—0.5, 0.5] x [—0.5, 0.5]. Black trajectories show island 
chains associated with various elliptic periodic orbits. There are no rotational invariant circles for this value of k. The red trajectory, 
started near the hyperbolic fixed point at (0.5, 0.0), is iterated 10 6 times and appears to be confined by rotational circles; however, it 
will escape after 2.3 x 10 9 steps. The blue and yellow trajectories are similarly only temporarily confined. This figure was generated 
with the StdMap application |Mei08| . 


2 It is convenient to choose this orientation to give conventionally oriented boundary integrals, 
orientation is positive. 


Note that that this means clockwise 
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Area-preserving maps are a special case of symplectic maps on a 2d-dimensional manifold that preserve a closed 
two-form. Symplectic maps are of interest because they arise as canonical Poincare maps for Hamiltonian flows 
|Mei92j . Thus the theory of (MMP84b] applies to the calculation of transport in a two degree-of-freedom Hamil¬ 
tonian system and has applications to mechanics, chemical reactions PP7!), IWWIOl IMS14] , plasma confinement 
[HM031 IBoo04| , celestial mechanics |WH911 IKLMROOl IMH01] , particle accelerator design [WilOlj , and condensed 
matter [T.TDF85] . More generally, volume-preserving maps arise as maps between sections or as stroboscopic maps 
for incompressible flows. Most physical applications correspond to three-dimensional flows such as the motion of a 
passive scalar in a fluid ICFP99 , IWQ041 IHall5l or mixing in granular media jSMQW08] . There are also nontrivial 
examples of periodically time-dependent two-dimensional flows, such as Aref’s blinking vortex |Are84j for which 
the area-preserving theory is directly applicable. 


A. Periodic Orbits and Resonance Zones 


Periodic orbits for area-preserving maps are typically born in pairs, one a saddle (hyperbolic) and the other 
elliptic. For example, the standard map 0 for k = 0 becomes the trivial shear (x,y) i-> (x + y,y); it integrable 
with action variable y. For this map, each circle on which y is rational is composed of period-?r orbits since 

( x n ,y„ ) = f n (%o, J/o) 

= (zo + ny 0 mod 1, y 0 ) = (x 0 , yo), 

when yo = — with m G Z and n G N. The minimal period of the orbit is n if gcd(m, n) = 1. These periodic orbits 
then have rotation number u = —, and we call them (to, n)-orbits. 

When the integrable shear is perturbed, at least two orbits survive for each rational rotation number; these are 
called the Birkhoff orbits, since they are predicted by the famous Poincare-Birkhoff theorem |Mei92l . In Fig. [I] such 
orbits can be seen at the centers of island chains. For the case of a twist map, like Q, Aubry-Mather theory shows 
that these orbits always exist; they are the so-called minimizing and minimax orbits because they can be obtained 
from a variational principle based on the action |Mat821 IAub83l . The minimizing orbit is hyperbolic whenever it is 
a nondegenerate minimum of the action, and the minimax orbit is elliptic when the perturbation is small enough 
IMM83I . 

An elliptic-hyperbolic pair of (to, ?r)-orbits forms the skeleton of a resonance zone consisting of n islands surround¬ 
ing the elliptic points with boundaries formed from the stable and unstable manifolds of the saddle IMMP871 IEas91| . 
One can view an island chain as a region of phase space in which the local rotation number is, roughly speaking, 
—. Computational methods for delineating resonances include Laskar’s frequency map [Las93j and the width of an 
orbit relative to a given rotation number |EMC93] . 

Geometrically, one can define a region of phase space that is bounded by broken separatrices of the saddle, an 
example is shown in Fig. [2j The boundary is formed using points on primary homoclinic orbits. Recall that a 
homoclinic point, 


to G W u {h) n W s {h), 


is primary if there are segments of stable, W s {h), and unstable, W u (h), manifold from the saddle h to to (that 
is, initial segments) that intersect only at their endpoints [Eas8611R.KW90] (see |LM00] for a generalization to 
higher dimensions). The resonance boundary in Fig. [2] consists of two initial segments of unstable manifold -for 
the “upper” and “lower” boundary—connected to two initial segments of stable manifold at primary homoclinic 
points. The region enclosed by the pair of broken separatrices is the resonance zone. The switch from unstable to 
stable manifold can occur at any point on the homoclinic orbit {toj}; even though the resulting geometry depends 
on this choice, the area of the resonance does not. In order for an orbit to leave the resonance zone it must pass 
through either the upper, E u or lower, £), exit set (we will say more about these in ) III A). Similarly, an orbit 
can enter the zone only by first landing in the preimage of one of the incoming sets, f~ i (I u ) or This 

construction can be done for any pair of (to, n)-orbits; a (1,3) resonance zone for Q is shown in Fig. [3] 

Each exit and incoming set is a lobe for the resonance, and the set E u U is a turnstile: it acts like a 

rotating door. 3 For each iteration of the map, trajectories escape through E u and arrive from / -1 (/„). 


! The name “turnstile” was suggested to us by Carl Murray. 
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FIG. 2: Resonance zone (shaded region) for the (0,1) orbits of 0 for k = 1.5. The boundary of the zone consists of two segments 
of unstable (red) and two segments of stable (blue) manifold of the (0,1) hyperbolic point joined at the primary homoclinic points 
mf. The upper and lower exit sets (dark green) and their first image (gray hatching) and preimage (light green) are shown, as are the 
preimage of the incoming sets (violet). 


B. Cantori 

The existence of invariant tori with sufficiently irrational frequency vectors in Hamiltonian systems that are 
sufficiently close to integrable is guaranteed by KAM theory |P5s01 . However, KAM theory says nothing about 
what happens to these structures far from integrability. The existence of remnants of these tori after they have 
been destroyed was first deduced by Serge Aubry and Ian Percival using variational arguments IAub781 lPer79l . 
Percival named the remnants cantori because—for an area-preserving map—they are invariant Cantor sets. Every 
orbit on a cantorus has the same irrational rotation number. 

Figure 0 shows the invariant circle of (|4|) with golden mean rotation number, u> = (j>, at the threshold of its 
destruction, k = k cr {(j)) ~ 0.97163540631 |Gre79| . and at seven larger values of k where it has become a cantorus. 
For twist maps, the projection of a cantorus onto the angle coordinate, x, is a Cantor set on the circle, so it has 
a countably infinite set of gaps. The orbit of each gap is a family of related gaps. For the standard map each 
cantorus appears to have one family, and each gap grows continuously from zero as k increases from k cr (oj). More 
generally, there may be more than one family of gaps and invariant circles may reform (BM93| . 

Mather and Aubry showed that for twist maps there is an invariant set for each irrational rotation number that 
is either an invariant circle or a cantorus IMat821 IAub83] . The same analysis applies to optical Hamiltonian flows 
[GoIOll . It is not clear, however, whether maps that do not satisfy a twist condition have cantori. Cantori do exist 
for a large class of maps close enough to an anti-integrable limit [AA901 !Aub95| . 

Questions I (Nontwist Cantori). What is the fate of invariant tori for area-preserving maps that do not satisfy a 
twist condition? Are there “shearless” cantori for nontwist maps like those studied by iFWAMOtfj ? 


The anti-integrable limit also can be used to show that certain higher-dimensional symplectic maps |MM92| and 
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FIG. 3: Resonance zone for the (1, 3)-orbits of the standard map when k = 1.5. The elliptic orbit (solid circles) is labeled (eo, ei, e 2 }, 
the hyperbolic (solid squares) {ho,hi,h, 2 }. The o symbols label a pair of (minimizing) homoclinic orbits (mj [ on the unstable (red) 
and stable (blue) manifolds of the saddle. The switch from unstable to stable manifold island occurs at the points mj and their two 
preimages. The shaded regions are turnstiles composed of an exit set, E u ( E {) (dark grey) and the preimage of an incoming set I u (/;) 
(light grey). 


higher-order recurrence relations modeling crystal lattices have analogues of cantori IKMR14I . However, Aubry- 
Matlrer theory has not been generalized to these cases. 

Questions II (Higher-Dimensional Cantori). When an invariant torus in a multi-dimensional symplectic map is 
destroyed, does it immediately become a cantorus, or are there other topological configurations such as Sierpinski 
curves ? The mechanism for destruction of two-tori in volume-preserving maps seems analogous to the symplectic 
case iFMlSty . but are there remnant cantori? 


III. THE PROBLEM OF TRANSPORT 

The inherent loss of predictability in chaotic dynamical systems implies that computation of individual trajecto¬ 
ries is not practical, nor is it especially useful. A theory of transport seeks to make predictions about the longtime 
evolution of ensembles of trajectories and answer questions such as “How long does it take, on average, to get from 
one region to another?”, or “What is the probability distribution of escape times from a region?” Thus for the 
chemical reaction AB + C —> A + BC , the first region might correspond to the separated reactants AB and C, 
and the second to the products A and BC. In plasma physics the first region may correspond to the interior of a 
confinement device such, as a tokamak, and escape to the divertor. 

The first step in a theory of transport is the computation of the flux across codinrension-one surfaces. For chaotic 
maps, this idea was used by Channon and Lebowitz |CL80I who computed escape times for Henon’s quadratic map. 
They used a broken separatrix of a period-five orbit to construct a simple, closed curve through which orbits can 
escape. In this section, I review the notions of flux and net flux for volume-preserving and exact volume-preserving 
maps. 

It is important to note that one of the fundamental measures of chaos, the Lyapunov exponent is not usually a 
good measure for transport. In |MMP84b] we said: 
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FIG. 4: Destruction of the golden mean invariant circle and creation of a cantorus for the standard map. At k = k cr (<f>) there is 
an invariant circle (top curve), but it is not smooth. Successive cantori, for k = 0.98 up to 1.04, are displaced downward for clarity. 
The largest gap forms around x = 0. The displayed orbits are actually saddle periodic orbits with rational rotation numbers that are 
convergents to the golden mean, see j)IV C| 


The Lyapunov exponents give the rate at which a nice cat in a region of phase space turns into a 
mixed-up cat. The turnstile rate constants give the mean rate at which bits of the cat are transported 
to regions of the phase space... 


A. Exit Sets and Flux 

A first step in a global theory of transport is to quantify the motion of an ensemble for short times by computing 
the flux of trajectories across codimension-one surfaces in phase space. Suppose R is a region with a piecewise 
smooth boundary dR. When R is invariant, f(R) = R, and dR 1 is a “complete” barrier. When R is “almost 
invariant” its boundary is a partial barrier: though trajectories may escape, the leakage is slow. For a flow, the flux 
through a surface is computed by integrating the normal component of the velocity over the surface IMM861 IMac94| . 
For a map, the flux is the volume that escapes from the region upon each iterate, i.e., the volume of the exit set 
E C R —the subset whose first image is not in R: 

E = {zGR | f(z) <£R} = {zeR\z(jt f~\R)} 

= R\f~\R ), 


( 6 ) 










see Fig. [5ja). The flux of escaping orbits from R is just the volume of the exit set, 


$(i?) = [i(E) = [ n. 

Je 

Similarly, the incoming set I for a region R is the set of points that have just entered 


( 7 ) 


I = R\f(R). 

Formally, R is almost invariant in the sense of measure if $(!?) -C One technique to compute almost 

invariant sets uses the Perron-Frobenius operator; this idea was pioneered by Dellnitz and Junge I)■)!)!) :. In lieu 
of this method, we will construct regions using resonance zones and cantori, for which the flux can be localized to 
turnstiles, see §IV| 



FIG. 5: Exit and Incoming Sets, (a) A region R , its exit set E with boundary U — S, and incoming set I. Orientations of the boundaries 
are shown by the arrows, (b) A rotational circle C on the cylinder and its image (dashed) and preimage (dotted). The net flux vanishes 
when the region above C and below /(C), labeled f(E) has the same area as that below C and above /(C), labeled I. 


A volume form ([2]) II is exact if 

f 1 = da (8) 

for some ( n — l)-form a. For example, the standard map Q is volume preserving with Q = dy A dx. Since x £ S 1 
is periodic, we cannot integrate f l with respect to this variable, but the Liouville form a = ydx, is well-defined on 
the cylinder and satisfies ([8|. When Q holds, Stokes’s theorem can be used to reduce volume integrals to surface 
integrals. 

The boundary of the exit set, dE, is the union of two oriented subsets U C dR and S C f~ 1 (dR) as sketched 
in Fig.^a). For a resonance zone, these are pieces of unstable and stable manifolds, recall Fig. [2] Assuming that 
the boundaries are chosen with an orientation consistent with that of R, then dE =U — S, and we can use ([8]) to 
rewrite 0 as 

^>(j R) = fi(E) = Ida = I a = fa — j a . (9) 

J e J dE Ju Js 

For two-dimensional maps, U and S are curves, and the line integrals in 0 are more efficient to compute than the 
original two-dimensional integral over E. 

In a sense, the complete transport problem is solved if one knows the orbits of E and / for each region of interest. 
For example, the minimum transit time from A to B is the first time that /*( Ea ) f~l Ib / 0. Similarly, a complete 
description of transport through a region requires knowing only the future history of the entering trajectories, i.e., 
the sets /*(/)■ These ideas form the basis for the theory of lobe dynamics as formulated by Rom-Kedar and Wiggins 
[RKW90IIRK94 . As we will recall in ^V| this is especially useful in the study of transport through resonance zones. 
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B. Exact Maps and Net Flux 


When a region R has finite measure, it is easy to see that the net flux through dR vanishes, so that the measure 
of the incoming and exit sets must be the same: 


m(J) = KR) - n f(R)) = y(R) - Kf~\R) n R) 
= v.{R)-n{Rnf-\R)) = n{E), 


where we used ([I]) and §• 

However, when the phase space is not contractible there may exist closed surfaces that that do not bound a 
region. For example the standard map 0 models a kicked rotor for which the canonical coordinate is an angle, 
and the phase space is the cylinder S x R. A rotational circle is a simple closed curve that is not contractible, 
for example Cq = {(x, 0) : x € §}. An area-preserving map on the cylinder need not have zero net flux through 
rotational circles; a simple example is the map f(x, y) = (x, y + 1) that shifts each rotational circle upward but still 
has det Df = 1. 

Zero net flux is equivalent to exactness: a map / is exact volume-preserving when there exists a Lagrangian form 
A such that 


f*a — a = dX. 

where a satisfies ([8]) (recall the discussion of pullback on page [2]). The Lagrangian for a two-dimensional 
a zero-form—a function —and is the discrete-time analog of the continuous-time Lagrangian of mechanics, 
for a flow, A is the time-integral of the Lagrangian along a trajectory between sections 1MM861 [WM97] , 

For example using a = ydx , the standard map 0 is exact with the Lagrangian 

A(x, y) = \y' 2 + jfi cos(27ra:), 

where y'{x,y) = y — sin(27r:r). This follows because 

f*a — a = y'dx' — ydx = y'd{x + y') — ydx 
= y'dy' + (y - ^ sin(27r x))dx - ydx 
= d{\y' 2 + ^2 cos(2ttx)), 

the differential of ©• In our original work for area-preserving twist maps we wrote A = F(x,x r ), which is one 
of the standard generating functions for canonical transformations IMMP84b] , I first learned of the more general 
formulation from Bob Easton jEas91j . 

For n-dinrensional maps A is an (n — 2)-form. For example, the ABC map [FKP88] . 

( x'\ f x + A sin 2 + C cos y \ 

y' J = I y + Lisina"' + A cos 2 | , (12) 

z' J \z + C sin y' + B cos x' J 


( 10 ) 

map is 
Indeed 

(ii) 


preserves the standard volume fl = dx A dy A dz, a fact that is most easily seen by viewing it as the composition 
of three shears, f = s z o s y o s Xl where, s x (x, y , z) = (x + A sin z + C cos y, y, z ), etc. Though it is natural to view 
(12) as a map on the three-torus, it can also be thought of as a map on the cover M = T 2 x 1. Now the volume 


form is exact with a = zdx A dy, and (12) is also exact with Lagrangian one-form 


A =A [sin(2)(d.T' + 
+ B sin (a ~!)dy' 


zdy) + cos {z)(dy — zdx')\ 
+ C cos(y')dx'. 


Here the primed variables are to be thought of as the functions defined in (12). Such a representation is especially 


pertinent when B,C <C 1, since then z is an action-like variable (invariant when B = C = 0), and one can use it 
to compute fluxes that correspond to vertical motion. 

It is also possible to generalize the idea of canonical generating functions to generating forms for volume¬ 
preserving transformations when a generalized twist condition is satisfied (LM09a| . 
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Questions III (Generating Forms). Can Lagrangian forms be used for perturbation theory and near-integrable 
normal forms, just as canonical generators are used for the symplectic case? Can they be used to build volume¬ 
preserving integrators, like those found for area-preserving maps in \ WM97( ? 

When / is homotopic to the identity, the net-flux “across” a closed (n — l)-dimensional surface S has meaning: 
it is the signed volume contained in the region R bounded by S and f(S). If dR = f(S) — S, Stokes’s theorem 
implies that 


$Net(S) = f n = f 

Jr Jf( 

= [ (fa - a) = [ dX = 0, 

Js Js 


lf(S)-S 


(13) 


because the integral of a total differential (exact form) over a closed surface S vanishes. For example in Fig. [5](b) , 
the rotational circle C on the cylinder is the upper boundary of an unbounded set with y —> — oo. This set has 
exit set E, the region below C whose image is above C, and incoming set J, the region below C whose preimage 
is above C. The region contained between f{C) and C is R = f{E) U I, and when the map has zero net flux (13) 
implies that p(f(E)) = p(E) = y(I). 

Maps with nonzero net flux also have interesting transport properties |FdlL15j . 


IV. ACTION AND FLUX 

Perhaps the most surprising aspect of the calculation of |MMP84b] is the connection between Mather’s action 
difference and the flux of trajectories across curves in phase space. 



FIG. 6: A three-dimensional lobe bounded by a pair of surfaces U (dark red) and S (light blue) that intersect along a loop 13 on their 
common boundaries. 


In the general context of an exact volume-preserving map, (101, let S and U be a pair of codimension one surfaces 
that together bound a region E, as sketched in Fig. [6j Suppose that 

B=UCS = dU = dS 


is the common (n — 2)-dimensional boundary of the surfaces. The volume of E is given by the difference between 
integrals over the surfaces as in ([9]). Each of these (n — l)-dimensional integrals can be partially integrated using 
(10) with Stokes’s theorem; for example, 


f a= f (—d\ + f*a) = — f A + f , 

Js Js Jb Jf(S) 


(14) 
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This iterative formula can be used to compute lobe volumes in a number of different contexts. For example, suppose 
that U = /"(£>); then we can iterate (14) n-times to obtain 


a 


= -£ 


A + 


a. 


t—0 


P(B) 


>f n (S ) 


(15) 


Note that the final integral in (15) is, by assumption, exactly the integral over the surface U, and so from we 
have 


p p n— 1 p 

H(E) = a- a = Y, I* 

JIA JS +—c\ J& 


A. 


(16) 


The last sum represents a generalized action of the periodic boundary B = f n (B). 

A second case in which (14) is useful is when S is a patch of stable manifold of some lower-dimensional invariant 


set, i.e., S C W S (A ), where A is a normally hyperbolic invariant manifold. In this case the limit of the remainder 
integral in (15) as n —» oo is zero, and thus 


a 


= -£ 




(17) 


is the forward action sum. Similarly, if U C W U (B) is a patch of unstable manifold of an invariant set B , it is 
appropriate to reverse (14) to iterate backwards, and obtain 


a 


= £ 


A. 


IfHB) 


(18) 


Lomeli and I showed how to compute flux for a 3D volume-preserving map |LM09b| . and MacKay showed 
how to compute flux over a saddle in a Hamiltonian flow }Mac89l IMS14] (a result implicit in |TJDF85] . as well). 
Nevertheless, the general problem in higher dimensions still has many mysteries. 

Questions IV (Multidimensional Flux). Can one use the flux formulas to compute flux through 2D manifolds in 
a 3D map, see e.g., IMJlAj . and through resonance zones formed from 2D normally hyperbolic invariant sets in a 
4 D symplectic map , see e.g., lGE91\j ? 


A. Flux Through Periodic Orbits 


When n = 2, B in (14) degenerates to a set of points, and f B A reduces to evaluation at its upper and lower 


limits. For example, suppose that U and S are curve segments joined at their endpoints, points on a pair of period-n 
orbits. An example would be e.g., the pair {eo} and {ho} that form an island chain, recall Fig. [3j In this case A is 
a function, and the boundary B = eo — ho is a pair of points, eo on the minimax orbit and h 0 on the hyperbolic 
orbit. Then (16) becomes 


n— 1 


KE) = £ [A(et) - A (h t )} = AW[e, h), 
(Periodic Orbit Flux), 


£=0 


(19) 


the difference between the actions of the two orbits. 

Indeed, AW[e, h\ is the flux “through” the pair of orbits, as can be seen by constructing (if possible) a region R 
with boundary C that satisfies two hypotheses. 

(HI). Suppose C = dR is a simple closed curve alternately joining successive points on a pair of period-n orbits, 
{e t } and {ht}. 

(H2). Suppose that /(C) DC = {e t } U {h t }. 
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An example of such a curve is sketched in Fig. [7] The requirement (H2) means that the exit and incoming lobes 
are bounded by segments of C and f(C) that connect pairs of points on the two orbits. By definition ([6]), f(E) is 
the region outside C and inside /(C). Each of the n-lobes of the image of the exit set is bounded by U t — S t where 
St (Ut) is the segment of C (f(C)) from h t to e t - Since f(S t ) = U t + 1 , and /(<S n _i) = Uq, (10) implies that the 
escaping flux is 


a 


n— 1 n, n— 1 p 

*(R) = ME)) = E / “ = E / f* a - 1 

t=0 JUt-St t=0 Js t 

n— 1 p n —1 

= E/ d\ = J2^- X ^)) = AW ^ h ^ 


—n JSt 


t—0 


t =0 


Thus the flux crossing C is exactly the action difference of (19). It is especially interesting, I think, that the flux 


is independent of the choice of curve C, given that it satisfies (HI) and (H2). Such curves are easy to construct 
for twist map examples, where the (to, n ) orbits are ordered relative to x , but it is not obvious that it is always 
possible to construct such curves for a pair of period n orbits. 



FIG. 7: Sketch of the construction of a partial barrier through a pair of (1, 3) orbits. The exited area is the region outside C and inside 
/(C), the dark (green) shaded areas. Under (HI) and (H2), this flux is (|19[l. 


Questions V (Periodic Flux). 


Under what conditions is it possible to find a curve C that satisfies (HI) and (H2)? 


B. Separatrix Lobe and Resonance Areas 


The iterative formula (14) can also be used to compute lobe and resonance zone areas for 2D maps. Recall from 
§11 A| that a resonance zone is bounded by segments of stable and unstable manifold of a hyperbolic periodic orbit 
joined at a primary homoclinic point. A separatrix lobe is a region bounded by segments of stable and unstable 
manifold that intersect on a neighboring pair of primary homoclinic points. 

More generally, following | Eas9lj . suppose that R is a region bounded by K alternating segments of stable and 
unstable manifold, thus dR = J2f=i(^ j + ^) where the lengths of /*(<S J ) and f^ l (W) limit to zero as t —> oo 
since they are stable and unstable segments, respectively. Denote the 2 K “corner points”, where these curves are 
joined, by to j and b? , j = 1,2,... K , so that dW = to/ — hf and dS J = /P +1 _ rn i ; with h K+1 = ft 1 . Labeling the 
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orbits of the corner points with subscripts, /‘(m- 7 ) = m J t , as usual, then (17) and (18) become 

r. OO 

/ . a = - Y ( A (^ +1 ) - A ( m t)) : 

J SI t=0 

f a = Y (a (mi) - A (hi)) . 

Xuj t—-oo 

Finally, summing over the K alternating segments to obtain the full boundary of R implies, from |9]), 


K 


K R ) = / «= Y 


t dR 


t =—oo j—1 


( 20 ) 


For example suppose that / and p are saddle fixed points, and that U C W u (f) and S C W s {p) are segments 
that intersect on a pair of heteroclinic points ho and mo as shown in Fig. [8j Then (20) applies with K = 1: there 
are just two segments and two corner points. Thus the area of the lobe is the action difference 


OO 

1*{R) = ^2 (A(m t ) - X(h t )) = AW[m, h], 

t — — OO 


( 21 ) 


(Heteroclinic Lobe Area). 

This applies to the computation of the lobe area shown in Fig. [2] The same formula applies when the saddles are 
periodic orbits. 


/'i 



I I 


FIG. 8: Lobe formed from segments of stable, S and unstable, IA , manifold of hyperbolic fixed points p and /, respectively, with 
heteroclinic orbits {fit} and {mt}. The lobe area p>(R) is the difference in between the actions of heteroclinic orbits, ©■ 


The action principle can also be used to compute the area of a resonance zone as well, recall §11 A[ Each island 
in the resonance is bounded by four segments of manifold and so the result ( |20| ) applies with K = 2. There are 
four corner points: h l and h r , the “left” and “right” saddle points of one island, and m , the “upper” and “lower” 
homoclinic points. For example, for the period-three resonance of Fig. [ 3 J the leftmost island has corners h l = / 12 , 
m 2 , h r = ho, and m^". Since each of the islands have the same area, the total area of the resonance is n times that 
of one island: 


OO 

p(R) = n (A(m+) - M h t) + x ( m 7 ) - X (K )), 

t =—OO 

(Period-n Resonance Area). 


( 22 ) 
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The result was first used to compute resonan ce areas in [MMP87] . Explicit results were obtained for the sawtooth 
map, a piecewise linear version of @ |CM8mlCDM+9f)j . 

For a twist map like ©> the rotational periodic orbits are vertically ordered according to their rotation numbers, 
and there is a resonance zone for each (to, n)-minimizing/minimax pair of periodic orbits. A choice of homoclinic 
points can be made so that different zones intersect only in their turnstiles, and since the net flux vanishes, the 
areas of the two intersections must be equal. Thus one can partition phase space into rotational resonance zones 
|MMP87] . An example is shown in Fig. [9j Qi Chen showed that when there are no rotational invariant circles, the 
total area of the infinitely many resonance zones is the full area of phase space, that is, the partition is “complete” 
|Che87j . One way to visualize this is to plot the area below each rotational curve as a function of rotation number 
IMMP87] . Rotational curves can be formed from invariant circles, cantori, or broken separatrices. This area 
function jumps by the value (22) at every rational rotation number. The partition is a complete devil’s staircase 
when the total rise of the staircase is the sum of the jumps. 


Questions VI (Resonance Partition). Is the resonance partition complete for nontwist area-preserving maps when 
there are no rotational invariant circles? Is there a similar (complete) partition for one-action, volume-preserving 
maps? How could one const ruct a res onance partition for AD symplectic maps (some of the intriguing geometry 
has been recently explored in jLRO + I^j)? 



FIG. 9: Partition of the phase space of the standard map into resonance zones for k = 1.1. Shown are portions of the stable (blue) 
and unstable (red) manifolds for the y, and y resonance zones. Exit sets (green) and Incoming sets (violet) for the y and y 

resonances shown. The partition is complete because there are no rotational invariant circles. 
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C. Cantori as Partial Barriers 


During the early years of the study of chaos, it was commonly observed that transport in a mixed phase space is 
highly nonuniform and is particularly slow as an orbit wends its way through regions containing recently destroyed 
tori. Indeed, Chirikov’s famous numerical experiments for the standard map seemed to show that the time to cross 
from y « 0 to y ~ 0.5 goes to infinity as k J, 0.989 )Chi79) : recall Fig. [l] Thus it seemed reasonable to speculate 
that cantori, even though they are typically hyperbolic, could act as strong barriers to transport. The importance 
of cantori was confirmed in | |MMP84b] when we discovered that the flux is locally minimal on cantori. We also 
were able to explain Chirikov’s result, showing that the average crossing time for © has a power law singularity 
with exponent 3.05 at k cr ((f> ) ss 0.971. 

For symplectic twist maps, every orbit on an invariant Lagrangian torus is a minimum of the action [MMS89] . 
For the area-preserving case, as noted in §11 B| a destroyed invariant circle becomes a cantorus: the minimizing 
orbits persist and are typically hyperbolic. Suppose that a cantorus T u with rotation number w has is a single 
family of gaps: there are two minimizing orbits {It}, {ry} C T u that correspond to the left and right endpoints of 
the gaps. When the cantorus is hyperbolic there exist c, A > 0 such that ||r t — Z t || < ce~ x ^. Aubry-Mather theory 
implies that there is another distinguished orbit {m t }, the minimax orbit, that falls in the gaps of the cantorus 
and is homoclinic to the orbit of gap endpoints, i.e., mt —» It as t —► ±oo. 

The flux through I u can be computed using the general action principle (201. Suppose C is any curve that 


satisfies hypotheses (HI) and (H2): it is a simple closed curve through the cantorus that crosses each gap going 
through the minimax point, and /(C) D C = T u U {m t } is exactly the cantorus and the minimax orbit. The same 
analysis as before implies that the flux crossing C is the action difference 


*(T U )= Y, (AW-A(l t ))=AW[m,l], 

t——oo 

(Cantorus Flux). 


(23) 


This action difference is also prominent in Aubry-Mather theory: when it is zero T u is an invariant circle, and when 
it is positive, T u is a cantorus |Mat82l !Aub83 :. 

Computing invariant tori is difficult since orbits with irrational rotation numbers are infinitely long. One algo¬ 
rithm uses a Newton-like method to compute the Fourier series for the conjugacy to rigid rotation of an invariant 
torus }HdlLS12l . However, this method does not seem to permit the computation of cantori. Another algorithm 
approximates a torus by periodic orbits. For the 2D case, this can be done using the continued fraction expansion 
w = [ao, ai,...], with ag € Z, aj £ N, of the rotation number to give a sequence of rationals 


n, 


— [ a 0, Ol) ^2, • • • fflj] 

= Oo + l/(Oi + 1/(02 + . . . + 1/Oj)), 


(24) 


the convergents of u>. A similar technique uses the Farey tree, which, beginning with an interval [//)-, bounded 
by neighboring rationals, m r ni — min r = 1 , constructs the daughter md = mi + m r , and nd = ni + n r , with 
F 1 < 7T 1 < 'n‘ L ' R- e P ea ti n g this construction on the two new intervals generates a binary tree that contains every 
rational in the original interval, and for which every irrational is a limit of an infinite path. 



1 1 2 3 5 8 13 

2’ 3’ 5’ 8’ 13’ 21’ 34’ "" 


When the Farey path limits to the rotation number of a cantorus, the periodic flux (19) limits to (23). Thus the 
strongest barriers, corresponding to the smallest flux, are associated with cantori. 

An alternative notion of minimal flux through cantori was formulated by Polterovich |Po w- 
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FIG. 10: Flux for the standard map at k = k,- r (/P) through periodic orbits on the Farey Tree. Reproduced with permission from 
Physica D 20:387 (1986) Copyright 1986 Elsevier Publishing. 


Questions VII (Cantori and Minimal Flux). Find a category of twist maps for which the flux is provably a 
monotone decreasing function on the Farey tree. For example, if the map is analytic, is the flux asymptotically 
monotone? 

Though there is no completely satisfactory generalization of continued fractions to higher dimensions, Kim and 
Ostlund invented a binary tree that has most of the properties of the single frequency Farey tree IK086I (it does 
not have a unique path to each rational vector). Even though there have been many numerical studies of the 
break-up of invariant tori for symplectic maps K.US9 . ITomQBt lCFL04| and volume-preserving maps IFM13] , there 
is no conclusive evidence for the formation of remnant tori that are Cantor sets. 

Questions VIII (Higher-Dimensional Cantori). Is there a multidimensional generalization of the Farey tree or 
continued fraction that is optimal for approximating tori and cantori? After destruction, is a torus replaced by a 
remnant invariant set with the same rotation vector; if so, what is its topology? 


V. EXIT AND RETURN TIMES 

Though a region’s escaping and entering fluxes give one measure of transport, fluxes represent only the first step, 
both dynamically and theoretically. More global measures are obtained from distributions of exit, crossing, and 
return times |RKW90i IMei97| . Early numerical experiments were done by Channon and Lebowitz, who showed 
that the exit time distribution from a resonance zone for the Henon map ([5]) has a surprising power law form [CL80I . 
This observation was greatly extended by Karney |Kar83| , who used integer arithmetic to eliminate floating point 
errors, and also by Chirikov and Shepelyansky |CS84I . From the many subsequent observations of power law decays, 
it seems clear that this phenomena is due to the presence of stable structures within the region of interest: these 
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structures are extremely sticky. However, there is still much controversy about whether the asymptotic decay is 
truly a power law, and— if it is- if there are universal values for the exponent. 


A. Exit Time Distributions 


Given a subset A of phase space M, the (forward) exit time of a point a £ A is the number of iterates for it to 
leave: 


t\(a) = min {t : f\a) £ A}. 


(25) 


Similarly, the backward exit time is 

*a(°) = min {t : /“‘(a) $. A}. 

In either case, the exit time is set to oo if the images remain in A for all t £ N. By definition § , the exit set, E, 
of a region is the set with forward exit time one, and the incoming set, /, is the set with backward exit time one. 
The transit time is the time for an orbit to move from the incoming to the exit set of A ; since this is constant along 
an orbit we can define 


t t r nsit (a)=t+(a)+r A (a)- 1 , 

subtracting the 1, to measure the time to get from I to E. Finally, the return time for a point a € A is 

t r f urn (a) = min{t : /‘(a) £ A}. 


(26) 

(27) 


Note that if a is not in the exit set, its return time is 1 since the minimum is taken over positive t. 

Chaos by its very definition precludes accurate long-time computations of (25)-(27). Nevertheless, one would like 


to obtain estimates of their distributions. Whenever A is has finite, nonzero measure, its exit time distribution , 


Pa O') = 


Ml a € A : *a( q ) =j}) 
MM 


(28) 


is the probability of exiting in time j. Another common measure of transport is the distribution of recurrence 
times: the probability that an orbit beginning in A has first return time j, 


return,,, A: Ef™n {a)= j}) 

Pa [J) ~ n(A) 


It is not hard to relate this distribution to others that we have defined. If a point in A is not in the exit set, then 
its return time is 1, so we have 


return ( i \ _ f*(A \E) ME) 

PA (1) - ~1W~ ~ ~ MAY 

Points that exit A are in E, and f(E) is the incoming set for M\A. The first return to A thus occurs when the 
orbit transits M \ A, so 


p r A etUrn (J) = 


MM) 

MM 


p"(i-1), 


j > i- 


(29) 


Thus the return time distribution to A is essentially the transit time distribution through M\A. 

To compute the statistics of exit and return times, any invariant subset of A should be excluded since it has 
infinite exit time. For example, the islands of stability that generically surround any elliptic periodic orbit should 
be excluded from exit time computations for the resonance zone of an area-preserving map. It is convenient to 
define the accessible subset of A is the subset of points that can be reached from the outside, or equivalently, the 
set of points with finite backward exit time: 


Aacc — {a £ A:t (a) < oo}. 
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In order that the distribution (28 1 be 


This set has the same volume as the set with finite forward exit time, 
normalized over N, it should be redefined using A acc instead of A. 

In 1947 Mark Kac proved a result for stochastic processes that can be translated into a result for deterministic 
maps |MMP84bl IMac941 IMei97| . Namely, suppose that M has finite measure and the invariant measure is normal¬ 
ized so that p(M) = 1. If M aC c is the subset of M that is accessible to orbits beginning in A, then the average first 
return time to A is 


/j.return\ _ 

Ua M = 


1 


li{A) 




% dp = 


n(M a 


Pi A ) 


Using the relation between transit time and return time distributions ( |29[ ), gives the equivalent result: 

Lemma 1 (Average Exit Time). If A is a measurable set with incoming set I and accessible set A a 
exit time for points in the incoming set is 


(30) 


the average 


(t+)i = = 


d( A acc 

dU) 


(31) 


One interesting aspect of Lem. T] is that it implies that the exit time distribution has finite mean, and by (29), that 
the mean return time to A is also finite when n(M) = 1. Higher moments of these distributions need not be finite; 
a counterexample is given by a simple shear map lMei97| . 

Since (29) is normalizable then one can define its cumulative distribution, the probability that the return is at 


least as long as k, 


P ireturn 
A 

This is often called the Poincare recurrence distribution (Zas85.. CS841ICS99I . The Poincare recurrence theorem 
states that when p(M) = 1, almost all trajectories return, i.e., Pff turn ( i) = l. 

In most cases there are no analytical formulae for exit and return time distributions. For a small perturbation to 
an integrable map with a homoclinic connection to a saddle, however, Melnikov theory can be used to estimate the 
exit set volume |RKLW90i |Wig92j . Melnikov theory, originally developed for flows, was formulated for maps by 
Easton |Eas84l lLMRR08j . Another case that can be treated is that of an adiabatic perturbation, where the lobes 
cover the region swept by the separatrices of the frozen time subsystems IKW91| . If the structures of the homoclinic 
tangle or trellis can be given a topological classification, like that developed by Easton IEas861 [Eas98| . then these 
methods can be extended to approximate additional intersection areas, and hence the exit time distributions 
fRK94l iMDOfil lMit09j . 

Questions IX (Higher-Dimensional Trellises). Can one classify trellises in 3 D volume-preserving and AD sym- 
plectic maps to obtain information about transport properties? 


(k) = Y J PA turn (j)- 

j=k 


B. Stickiness and Anomalous Diffusion 


If a map is uniformly hyperbolic then the number of periodic orbits of a given period grows exponentially at a 
rate given by the topological entropy. This results in an exponential form for exit time distributions [ZT911IHSV99I : 


PaU) 


o-j/t 


However, as was first observed by |CL80| . such distributions for area-preserving maps with mixed regular and 
chaotic orbits appear to decay much more slowly: 


PaU) ~ t 


- 7-1 


(32) 


Such algebraic decay is a signal of the stickiness of stable structures in the phase space, see §VI| 

The existence of the average (31) implies that 7 > 0. It can be shown that the exponent for / and / -1 are the 


same |Eas93] . thus (26) implies that the transit time distribution will have this exponent as well, and hence, by 
(29), so will the return distribution. Consequently, the cumulative recurrence distribution will decay as 




( 33 ) 
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Many studies have found 7 ~ 1.5 |Ven09j . though there are fluctuations that persist for large times [CA13j , and 
it is still controversial whether the power law (331 is valid asymptotically, and if it is, whether there are classes of 
dynamical systems for which the exponent 7 is universal. 


Questions X (Power Law Decay). Do exit and return-time distributions in typical volume-preserving maps with 
mixed regular and chaotic components have an asymptotic power law form? Is there a universal exponent 7 ? 


The exponent 7 is related to another measure of transport, the growth of the mean square displacement (||;r t — 
x 0 || 2 ) with time. For Brownian motion, the squared displacement grows linearly in time and there exists a diffusion 
coefficient, 


D = „ lim ^f(\\ x T-x 0 \\ 2 ). 

T—too 11 

However, it is observed that maps like @ sometime exhibits anomalous momentum diffusion with 

{(yt - y 0 ) 2 ) ~ 1». 


(34) 


Super-diffusion, i.e., (3 > 1, is observed in the standard map when there are accelerator modes, stable orbits that 
advance linearly in the momentum direction |Chi791 IKRW821IRKZ981 IVenflSl IMR.141 IAR.14] . An accelerator island 
will drag nearby chaotic orbits along, so that even when the ensemble average does not include the accelerating 
island, the effect is that j3 > 1. The simplest accelerator mode forms at k = 2i r, where after one step, (xi,yi) = 
(xo, 2 /o + l)i however, there are probably high-period accelerator modes arbitrarily close to k cr {(j>). More general 
maps that do not have the vertical periodicity of the standard map do not have accelerator modes; however, the 
effect of sticky islands is still important. 

Karney |Kar83| showed that (33) and (34) can be related through the force correlation function, C(t). Indeed 
the effective diffusion after T steps is 


^((yT-yo) 2 ) = kC(o) + J 2 c (t). 


t —1 


Karney’s result is that the power law (33) implies C{t) ~ t 7+1 , which then, by (34), gives 

7 + /3 = 3. 


(35) 


Consequently, the algebraic decay of the recurrence distribution and anomalous diffusion are two sides of one coin. 
Since numerical experiments estimate that 7 ss 1.5, then (35) implies that fj ss 1.5 as well. This has been confirmed 


by Venegeroles using the Perron-Frobenius operator to analytically estimate diffusion due to accelerator modes 
[VenOSl rVenf)9] . 

Stickiness can also be studied using finite time Lyapunov exponents: the distribution of exponents is bimodal 
due to orbits sticking near elliptic regions [SLV05] . This is also observed in higher-dimensional systems |MBR14j . 


VI. MARKOV MODEL FOR TRANSPORT 

The phase space of a typical area-preserving map is a complex mixture of islands and chaotic orbits, recall Fig. [I] 
Each island, a collection of nested invariant circles, is also surrounded by periodic orbits that have rational rotation 
numbers relative to the elliptic periodic orbit at the island’s center. Each elliptic periodic orbit also generically 
forms islands. Numerically, there appear to be infinitely many islands comprising a fat fractal: a chaotic region has 
nonzero area but its boundary has non-integer dimension |UF85| . Thus, embedded in the chaotic sea is an infinite 
hierarchy of islands around islands, see Fig. [TT] 

A connected chaotic component can be partitioned into tree of states (MMP84b| IMQ 86 | . The nodes of the tree 
correspond to regions bounded by partial barriers. Each state is a portion of the chaotic component surrounding 
a given elliptic island, i.e., the accessible subset of a resonance zone. Transport between states is mediated by 
turnstiles. These could be lobes in resonance boundaries; however, at least for twist maps, there are also cantori 
between resonances (there are irrationals between every pair of rationals). Since lowest flux is that through a 
cantorus, recall §IV C[ these form the most restrictive partial barriers. The elliptic orbits inside a resonance zone 
also form islands that are separated by cantori. 
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FIG. 11: Island around island structure for the Henon map § at a = —0.17197997940. Successive enlargements are indicated by the 
arrows. Boundary circles for the states S = 1 and 10 are shown (brown). For this parameter, the map has a self-similar sequence of 
islands around islands, each with relative period 7. Five levels of this sequence are shown. Reproduced with permission from Physical 
Review E 90(6):062923 (2014). Copyright 2014 AIP Publishing. 


Transport through a chaotic region appears to be anomalous whenever it has a boundary component that is a 
stable invariant set. For area-preserving maps these are typically boundary circles of islands, the outermost of the 
encircling loops. It is these that give rise to the power law (32). 
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The simplest model contains a single island so that the chaotic component is bounded by a single invariant 
circle. The stickiness of a boundary circle is caused by a sequence of cantori that limit on the boundary with 
ever-decreasing flux. As John Greene explained to us in 1983 |MMP84b] . such a scenario yields a power law (|32|; 
however, the exponent is much larger than that observed: 7 = 3.05 |HCM85| or 7 = 3 |CS99j . Alternatively, if only 
the hierarchy of islands around islands is kept, ignoring the cantori, the exponent is still too large: 7 = 2.25 |ZEN97| . 
The insufficiency of either model was also confirmed by the observations of IWHK031 : a trajectory following either 
one of these hierarchies has an increasing probability of becoming trapped by secondary islands or by smaller flux 
cantori as time increases. This phenomena results in transient local density accumulations jSS87llMei94| . Note also 
that systems not smooth enough to have cantori and with no hierarchy of islands, such as some piecewise linear 
maps, have a different exponent, 7 = 2 [AMKOfij 

The simplest partition that keeps both phenomena is a binary Cayley tree (nodes have a parent and two daugh¬ 
ters), see Fig. 12 Each node in the tree is labeled by a sequence S = sqSi ..., with € {0,1}. The root of the 


tree, denoted 0, can be taken to be an absorbing state : whenever a trajectory reaches this state, it is deemed to 
have escaped. For the Henon map of Fig. El the null state is the region outside the resonance zone formed by the 
stable and unstable manifolds of the hyperbolic fixed point (blue and red curves). 

Transitions to the right and left on the tree correspond to distinct topological structures in phase space. A 
transition to the right corresponds to motion towards a given island boundary. The generically irrational rotation 
number of a boundary circle has a continued fraction expansion, and each of its convergents, (24), corresponds to 
a periodic orbit |GMS86j . The resonance zones of every other convergent lie in the chaotic component, giving an 
infinite sequence of levels that limit on the boundary circle. A transition from one level to the next corresponds to 
adding a “1” to the state, S —> 51. 

Each elliptic island within a state 5 also has its own boundary circle and corresponding infinite sequence of levels 
of encircling periodic orbits. We denote the outermost of these new levels by 50; thus a transition to the left on 
the tree, S —>• 50, corresponds to being trapped around an island of higher class | MMP84bl IMei86| . For example, 
in Fig. 11 the outermost elliptic islands in the upper right pane form a (1,7) chain. This represents state 5 = 1 in 
Fig. 12 One island in this chain is enlarged in the upper left panel of Fig. [llj it is also encircled by a (1, 7) chain 
that has rotation number ^ relative to the seventh image, /', so that this orbit has period 49 relative to /. In the 
figures this island is designated by its rotation number sequence ^ ■ J- These are the outermost elliptic islands of 
class one corresponding to the state 5 = 10. 

If there is only one significant elliptic island in each state, then the tree is binary. More generally, some states in 
the tree can have more daughters when there are multiple islands in a given state, and some may have fewer if the 
periodic orbits are hyperbolic |AFM14| . 



FIG. 12: Tree for the Henon map § at a = —0.17197997940 with the phase space shown in Fig. m where it has a period-seven 
self-similar hierarchy, from [AFM14] . Each node is labeled by the state S and the winding number sequence for the main island of that 
class, 21 : P 2 ; Several transition probabilities, Ps^-S' > are a l so indicated. Reproduced with permission from Physical Review E 
90(6):062923 2 (2014). Copyright 2014 AIP Publishing. 
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The simplest approximate model for transport on the tree of states is Markov. The assumption is that the chaos 
within a state S causes rapid mixing, and so the probability of a trajectory leaving at any time is proportional to 
the flux, of trajectories from 5 to state S', where S' is the parent, DS, or a daughter, SO or 51. This 

flux is the area of the turnstile of the most resistant cantorus that divides the states and can be computed as an 
action difference, recall f IV Note that the flux is symmetric, A Ws s’ = A Ws',s, when the map is exact, recall 
This model seems reasonable when the average exit time (31) is large compared to a Lyapunov time—the 


t III B 


time for significant memory loss in the dynamics. However, for a near-integrable system, the mixing approximation 
probably will not work, and there are strong correlations that result in oscillatory exit time distributions |RK94| . 


Questions XI (Markov Transport Models). Is there a limit in which the Markov model is valid? When can 
higher-order correlations be neglected? 


Thus a Markov model for transport on the tree is defined by transition probabilities ps->S' for each pair of 
connected nodes on the tree, recall Fig. [12} Under the assumption that a trajectory has equal probability to be 
anywhere in 5, the probability of such a transition is 


PS—*S' 


AfUs.s' 

As 


(36) 


where A$ is the accessible area of 5. By Lem. |T| (|36} is exactly the inverse of the average transit time through the 
region. Following [MQ 86 j . it is convenient to categorize the change in transition probabilities on the tree by two 
ratios, 


i(S) = 


PS^Si 
PS^DS ’ 


£i(5) = 


PSi^S 
PS^DS ’ 


(37) 


for i £ {0,1}. The first ratio measures asymmetry between motion “down” and “up” the tree, and the second 
measures the slowing of time scales deeper in the tree. 

For example, renormalization theory implies that for a noble invariant circle, W\ = 0.053112: it is 19 times more 
likely for an orbit to move away from the circle—to the parent state—than to move one level deeper )Mac83| . For 
the self-similar sequence of period-seven islands shown in Fig. EL the class ratio is even smaller, W 2 = 0.014158, 
implying that trapping around an island is an even rarer event IMei 86 l . For this reason, anomalous transport is 
observable only for extremely long trajectories. Accordingly, it is debatable whether an asymptotic regime is reached 
even for the longest computations (e.g., quadruple precision computations for orbits of length 10 12 |WHK03j ). 

In the original Markov tree model, Ed Ott and I assumed that the tree is self-similar [MQ 86 I . This would 
asymptotically hold if every boundary circle were noble and every island chain had the same period relative to its 
parent. In this case the ratios (37) are independent of the state 5, though they depend on the choice of class, * = 0, 
or level, i = 1. More generally, the results of [ GMS 86 ] imply that the level hierarchy is approximately self-similar, 
and that for almost all boundary circles 


wi ~ £i 05 . 


(38) 


Similarly, there are cases in which the class hierarchy is asymptotically self-similar |Mei 86 l . Numerical studies again 
show that the scaling coefficients are related, 


wo - £ o 19 > 


(39) 


even as the period of the successive islands varies. 

If the binary tree is exactly self-similar, then there is an integral equation for the recurrence distribution (MQ861 . 
A solution to this gives the power law (32), with an exponent determined by the dispersion relation 


VJ 0 £ 0 J 


■ Wi£ 1 7 


= 1. 


Given the relations, (38)-(39), there are still two parameters needed to solve for 7 . Reasonable values for the e. 


give 7 = 1.96, which is closer to the observed value, 7 ~ 1.5, than the models with no tree hierarchy. The analysis 
of [MO 86 ] implies that each branch added to the tree decreases the value of 7 . 

Note, however, that for a typical map, neither self-similar scenario will hold 
depend upon the state 5. 


and the scaling factors ( [37} will 

Recent studies have modeled the lack of self-similarity statistically. 


The idea is that 

an orbit will sample different regions of phase space, and its long-time behavior will be equivalent to an ensemble 
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of scalings. Effectively, the coefficients of the Markov tree model should be taken as a random draw from some 
probability distribution. Numerical studies do appear yield different values of 7 in different regions of phase space 
IWHK031 IAdC13] . Using such a model, Cristadoro and Ketzmerick find 7 = 1.57, in remarkable agreement with 
numerical experiments [CK08] . A similar model, using a different ensemble was studied in [CA13] . Indeed, it is 
not clear which random ensemble of scaling coefficients is valid, though there may be some sense in which there is 
a universal ensemble IA F.Yl 1 l] . The ensemble picture also provides an explanation for the difficulty of observing 
an asymptotic regime: the correlation function of the fluctuations in the mean decay exponent have a log-periodic 
behavior that decays very slowly in time lCA13j . 

The Markov model, then, not only gives a qualitative explanation for the power law decay of the exit time 
distribution (321, but also gives a value of 7 close to that observed. Of course, this is only for the area-preserving 
case. 


Questions XII (Higher-Dimensional Transport). Are there Markov tree models that quantitatively explain trans¬ 
port in higher-dimensional volume-preserving and symplectic maps? 


VII. THE NEXT THIRTY YEARS 

I have tried to summarize both the state of the art in the theory of transport for conservative dynamical systems 
and to indicate some directions that seem to me to be still ripe for future research. 

There are other fruitful techniques that I did not include in this survey, such as the characterization of mixing 
in two-dimensional flows by the braiding of trajectories [BASOOi IFT11I . and the set theoretic methods based on 
the Perron-Frobenius or Koopman operator formalisms 1-Vo!)!) i! Fro05l IBMM121 IBS 13] . I also did not discuss the 
continuous time and space approximation of anomalous transport models to obtain phenomenological, fractional 
kinetic equations. This theory, which is reviewed in IZEN971 IZas02j . does not have as close a connection with the 
underlying topology of the dynamics as the discrete Markov tree. Finally, I did not discuss the application of the 
ideas of turnstiles to semi-classical and quantum mechanics. Indeed, in |MMP84b] we speculated that classical 
cantori whose turnstiles were small in units of Planck’s constant might play a significant role in the semiclassical 
mechanics of molecules, and MacKay and I later applied this idea to compute the ionization threshold for hyd rogen 
in a strong microwave field [MM 88 j . Several recent advances in this area confirm our original speculation M BK + 12l 
IRT BK1 4| . 

There is much that remains intriguing about this problem: the infinite complexity of phase space and the long¬ 
time correlations that result in apparent power law decays and anomalous diffusion. Even after thirty years, I still 
cannot resist the allure of phase space portraits for simple area-preserving maps. 

In [MMP84F we noted 

The theory is not yet complete, but we believe that the ideas introduced here are of central importance 
for a broad range of systems that are neither completely ordered nor completely chaotic. 

While there has been considerable progress, I think—as indicated by the series of questions posed in this paper—that 
there is still much to do in the study of transport. 


Acknowledgments 

The work was partially supported NSF grant DMS-1211350. I would like to thank Robert MacKay, Ian Percival 
and Bob Easton for long and fruitful collaborations, many discussions, and a profusion of insights. 


[AA90] S. Aubry and G. Abramovici. Chaotic trajectories in the standard map, the concept of anti-integrability. Physica 
D, 43:199-219, 1990. http://dx.doi.org/10.1016/0167-2789(90)90133-A 
[AdC13] C.V. Abud and R.E. de Carvalho. Multifractality, stickiness, and recurrence-time statistics. Phys. Rev. E, 
88(4):042922, 2013. http://link.aps.org/doi/10.1103/PhysRevE.88.042922 
[AFM14] O. Alus, S. Fishman, and J.D. Meiss. Statistics of the island-around-island hierarchy in Hamiltonian phase space. 
Phys. Rev. E, 90(6):062923, 2014. http://link.aps.org/doi/10.1103/PhysRevE.90.062923 




































24 


[AMK06] E.G. Altmann, A.E. Motter, and H. Kantz. Stickiness in Hamiltonian systems: From sharply divided to hierarchical 
phase space. Phys. Rev. E, 73(2):026207, 2006. http://dx.doi.org/10.1103/PhysRevE.73.026207 
[AR14] T. Albers and G. Radons. Weak ergodicity breaking and aging of chaotic transport in Hamiltonian systems. Phys. 

Rev. Lett., 113(18):184101, 2014. http://link.aps.org/doi/10.1103/PhysRevLett.113.184101 
[Are84] H. Aref. Stirring by chaotic advection. J. Fluid Mech., 143:1-21, 1984. http://dx.doi.org/10.1017/ 

S0022112084001233 

[ASZ91] V.V. Afanasiev, R.Z. Sagdeev, and G.M. Zaslavsky. Chaotic jets with multifractal space-time random walk. Chaos, 
1(1):143-159, 1991. http://dx.doi.Org/10.1063/l.165824 

[Aub78] S. Aubry. The new concept of transitions by breaking of analyticity in a crystallographic model. In A.R. Bishop 
and T. Schneider, editors, Solitons and Condensed Matter Physics, volume 8 of Solid State Sciences, pages 264-277. 
Springer-Verlag, New York, 1978. http://dx.doi.org/10.1007/978-3-642-81291-0 
[Aub83] S. Aubry. The twist map, the extended Frenkel-Kontorova model and the devil’s staircase. Physica D, 7:240-258, 

1983. http://dx.doi.org/10.1016/0167-2789(83)90129-X 

[Aub95] S. Aubry. Anti-integrability in dynamical and variational problems. Physica D, 86:284-296, 1995. http://dx.doi. 
org/10.1016/0167-2789(95)00109-H 

[BASOO] P.L. Boyland, H. Aref, and M.A. Stremler. Topological fluid mechanics of stirring. J. Fluid Mech., 403:277-304, 
2000. http://dx.doi.org/10.1017/S0022112099007107 

[BK84] D. Bensimon and L.P. Kadanoff. Extended chaos and disappearance of KAM trajectories. Physica D, 13:82-89, 

1984. http://dx.doi.org/10.1016/0167-2789 (84)90271-9 

[BM93] C. Baesens and R.S. MacKay. Cantori for multiharmonic maps. Physica D, 69(l-2):59-76, 1993. http://dx.doi. 
org/10.1016/0167-2789(93)90180-9 

[BMM12] M. Budisic, R. Mohr, and I. Mezic. Applied Koopmanism. Chaos, 22(4):047510, 2012. http://dx.doi.org/10. 
1063/1.4772195 

[Boo04] A.H. Boozer. Physics of magnetically confined plasmas. Rev. Mod. Phys., 76:1071- 1141, 2004. http://link.aps. 
org/doi/10.1103/RevModPhys.76.1071 

[BS13] E.M. Bollt and N. Santitissadeekorn. Applied and Computational Measurable Dynamics, volume 18 of Mathematical 
Modeling and Computation. SIAM, Philadelphia, 2013. 

[CA13] R. Ceder and O. Agam. Fluctuations in the relaxation dynamics of mixed chaotic systems. Phys. Rev. E, 
87(1):012918, 2013. http://dx.doi.org/10.1103/PhysRevE.87.012918 
[CDM + 90] Q. Chen, I. Dana, J.D. Meiss, N. Murray, and I.C. Percival. Resonances and transport in the sawtooth map. 

Physica D, 46:217-240, 1990. http://dx.doi . org/10.1016/0167-2789(90)90037-P 
[CFL04] A. Celletti, C. Falcolini, and U. Locatelli. On the break-down threshold of invariant tori in four dimensional maps. 

Regul. Chaotic Dyn., 9(3):227-253, 2004. http://dx.doi.org/10.1070/RD2004v009n03ABEH000278 
[CFP99] J.H.E. Cartwright, M. Feingold, and O. Piro. An introduction to chaotic advection. In H. Chate, E. Villermaux, 
and J.M. Chomez, editors, Mixing: Chaos and Turbulence, pages 307-342. Kluwer, 1999. 

[Che87] Q. Chen. Area as a devil’s staircase in twist maps. Phys. Lett. A, 123:444-450, 1987. http://dx.doi.org/10.1016/ 
0375-9601(87)90343-4 

[Chi79] B.V. Chirikov. A universal instability of many-dimensional oscillator systems. Phys. Rep., 52:265-379, 1979. http: 
//dx.doi.org/10.1016/0370-1573 (79)90023-1 

[CK08] G. Cristadoro and R. Ketzmerick. Universality of algebraic decays in Hamiltonian systems. Phys. Rev. Lett., 
100:184101, 2008. http://link.aps.org/doi/10.1103/PhysRevLett.100.184101 
[CL80] S.R. Channon and J.L. Lebowitz. Numerical experiments in stochasticity and heteroclinic oscillation. In R.H.G. 

Helleman, editor, Nonlinear Dynamics, volume 357, pages 108-118. Ann. New York Acad. Sci., New York, 1980. 
[CM89] Q. Chen and J.D. Meiss. Flux, resonances and the devil’s staircase for the sawtooth map. Nonlinearity, 2:347-356, 
1989. http://dx.doi.org/10.1088/0951-7715/2/2/009 

[CS84] B.V. Chirikov and D.L. Shepelyansky. Correlation properties of dynamical chaos in Hamiltonian systems. Physica 
D, 13:395-400, 1984. http://dx.doi.org/10.1016/0167-2789(84)90140-4 
[CS99] B.V. Chirikov and D.L. Shepelyansky. Asymptotic statistics of Poincare resurrences in Hamiltonian systems with 
divided phase space. Phys. Rev. Lett., 82:528-531, 1999. http://dx.doi.org/10.1103/PhysRevLett.82.528 
[DJ99] M. Dellnitz and O. Junge. On the approximation of complicated dynamical behaviour. SIAM J. Numer. Anal, 
36(2):491-515, 1999. http://dx.doi.org/10.1137/S0036142996313002 
[Eas84] R.W. Easton. Computing the dependence on a parameter of a family of unstable manifolds: Generalized Melnikov 
formulas. Non. Anal. Th. Meth. and AppL, 8:1-4, 1984. http://dx.doi.org/10.1016/0362-546X(84)90023-3 
[Eas86] R.W. Easton. Trellises formed by stable and unstable manifolds in the plane. Trans. Am. Math. Soc., 294(2):719-732, 
1986. http://www.jstor.org/stable/2000211 

[Eas91] R.W. Easton. Transport through chaos. Nonlinearity, 4:583-590, 1991. http://dx.doi.org/10.1088/0951-7715/ 
4/2/017 

[Eas93] R.W. Easton. Transport of phase space volume near isolated invariant sets. J. of Dyn. and Dijf. Eq., 5(3):529-535, 
1993. http://dx.doi.org/10.1007/BF01053534 

[Eas98] R.W. Easton. Geometric Methods for Discrete Dynamical Systems. Cambridge University Press, Cambridge, 1998. 
[EMC93] R.W. Easton, J.D. Meiss, and S. Carver. Exit times and transport for symplectic twist maps. Chaos, 3(2):153-165, 




25 


1993. http://dx.doi.org/10.1063/1.165981 

[FdlL15] A.M. Fox and R. de la Llave. Barriers to transport and mixing in volume-preserving maps with nonzero flux. 

Physica D, 295-296:1-10, 2015. http://dx.doi.org/10.1016/j .physd.2014.12.003 
[FKP88] M. Feingold, L.P. Kadanoff, and O. Piro. Passive scalars, three-dimensional volume-preserving maps and chaos. J. 

Stat. Phys., 50:529, 1988. http://dx.doi.org/10.1007/BF01026490 
[FM13] A.M. Fox and J.D. Meiss. Greene’s residue criterion for the breakup of invariant tori of volume-preserving maps. 

Physica D , 243(l):45-63, 2013. http://dx.doi.Org/10.1016/j.physd.2012.09.005 
[Fro99] G. Froyland. Using Ulam’s method to calculate entropy and other dynamical invariants. Nonlinearity, 12:79-101, 
1999. http://dx.doi.Org/10.1088/0951-7715/12/l/ 

[Fro05] G. Froyland. Statistically optimal almost-invariant sets. Physica D, 200(3-4):205-219, 2005. http://dx.doi.org/ 
10.1016/j.physd.2004.11.008 

[FT11] M.D. Finn and J-L. Thiffeault. Topological optimization of rod-stirring devices. SIAM Review, 53(4):723-743, 2011. 
http://dx.doi.org/10.1137/100791828 

[FWAM06] K. Fuchss, A. Wurm, A. Apte, and P. J. Morrison. Breakup of shearless meanders and “outer” tori in the standard 
nontwist map. Chaos, 16:033120, 2006. http://dx.doi.Org/10.1063/l.2338026 
[GE91] R.E. Gillilan and G.S. Ezra. Transport and turnstiles in multidimensional Hamiltonian mappings for unimolecular 
fragmentations: Application to van der Walls predissociation. J. Chem. Phys., 94(4):2648-2668, 1991. http://dx. 
doi.org/10.1063/1.459840 

[GMS86] J.M. Greene, R.S. MacKay, and J. Stark. Boundary circles for area-preserving maps. Physica D, 21:267-295, 1986. 
http://dx.doi.org/10.1016/0167-2789(86)90005-9 

[GolOl] C. Gole. Symplectic Twist Maps: Global Variational Techniques, volume 18 of Advanced Series in Nonlinear Dy¬ 
namics. World Scientific, 2001. 

[Gre79] J.M. Greene. A method for determining a stochastic transition. J. Math. Phys., 20:1183-1201, 1979. http: 
//dx.doi.org/10.1063/1.524170 

[Hall5] G. Haller. Langrangian coherent structures. Annual Review of Fluid Mechanics, 47:137-162, 2015. http://dx.doi. 
org/10.1146/annurev-fluid-010313-141322 

[HCM85] J.D. Hanson, J.R. Cary, and J.D. Meiss. Algebraic decay in self-similar Markov-chains. J. Stat. Phys., 39(3-4):327- 
345, 1985. http://dx.doi.org/10.1007/BF01018666 

[HdlLS12] G. Huguet, R. de la Llave, and Y. Sire. Computation of whiskered invariant tori and their associated manifolds: 

New fast algorithms. Disc. Cont. Dyn. Sys., 32(4):1309-1353, 2012. http://dx. doi. org/10.3934/dcds. 2012.32.1309 
[Hen69] M. Henon. Numerical study of quadratic area-preserving mappings. Q. J. Appl. Math., 27:291-312, 1969. 

[HM03] R.D. Hazeltine and J.D. Meiss. Plasma Confinement. Dover Publications, Mineola, NY, 2nd edition, 2003. 

[HSV99] M. Hirata, B. Saussol, and S. Vaienti. Statistics of return times: A general framework and new applications. Comm. 

Math. Phys., 206(l):33-55, 1999. http://dx.doi.org/10.1007/s002200050697 
[Kar83] C.F.F.K. Karney. Long time correlations in the stochastic regime. Physica D, 8:360-380, 1983. http://dx.doi. 
org/10.1016/0167-2789(83)90232-4 

[Kec67] J.C. Keck. Variational theory of reaction rates. In Advances in Chemical Physics, pages 85-121. John Wiley & 
Sons, Inc., 1967. http://dx.doi.org/10.1002/9780470140154.ch5 
[KLMROO] W.S. Koon, M.W. Lo, J.E. Marsden, and S.D. Ross. Heteroclinic connections between periodic orbits and 
resonance transitions in celestial mechanics. Chaos, 10(2):427-469, 2000. http://dx.doi.Org/10.1063/l.166509 
[KM89] H.T. Kook and J.D. Meiss. Periodic-orbits for reversible, symplectic mappings. Physica D, 35(l-2):65-86, 1989. 
http://dx.doi.org/10.1016/0167-2789(89)90096-1 

[KMR14] V. Knibbeler, B. Mramor, and B. Rink. The laminations of a crystal near an anti-continuum limit. Nonlinearity, 
27(5):927-952, 2014. http: //dx. doi . org/10.1088/0951-7715/27/5/927 
[KO 86 ] S. Kim and S. Ostlund. Simultaneous rational approximations in the study of dynamical systems. Phys. Rev. A, 
34:3426-3434, 1986. http://dx.doi.org/10.1103/PhysRevA.34.3426 
[KRW82] C.F.F.K. Karney, A.B. Rechester, and R.B. White. Effect of noise on the standard mapping. Physica D, 4:425-438, 
1982. http://dx.doi.org/10.1016/0167-2789(82)90045-8 

[KW91] T.J. Kaper and S. Wiggins. Lobe area in adiabatic Hamiltonian systems. Physica D, 51:205-212, 1991. http: 
//dx.doi.org//10.1016/0167-2789(91)90233-Y 

[Las93] J. Laskar. Frequency analysis for multi-dimensional systems. Global dynamics and diffusion. Physica D, 67:257-283, 
1993. http://dx.doi.org/10.1016/0167-2789(93)90210-R 

[LM00] H.E. Lomelf and J.D. Meiss. Heteroclinic primary intersections and codimension one Melnikov method for volume¬ 
preserving maps. Chaos, 10(1): 109—121, 2000. http://link.aip.Org/link/7CHA0EH/10/109/l 
[LM09a] H.E. Lomelf and J.D. Meiss. Generating forms for exact volume-preserving maps. Disc. Cont. Dyn. Sys. Series S, 
2(2):361-377, 2009. http://dx.doi.Org/10.3934/dcdss.2009.2.361 
[LM09b] H.E. Lomelf and J.D. Meiss. Resonance zones and lobe volumes for volume-preserving maps. Nonlinearity, 22:1761- 
1789, 2009. http://dx.doi.Org/10.1088/0951-7715/22/8/001 
[LMRR08] H.E. Lomelf, J.D. Meiss, and R. Ramfrez-Ros. Canonical Melnikov theory for diffeomorphisms. Nonlinearity, 
21:485-508, 2008. http://dx.doi.Org/10.1088/0951-7715/21/3/007 
[LRO + 14] S. Lange, M. Richter, F. Onken, A. Backer, and R. Ketzmerick. Global structure of regular tori in a generic 4D 



26 


symplectic map. Chaos, 24:024409, 2014. http://dx.doi.Org/10.1063/l.4882163 
[Mac83] R.S. Mac.Kay. A renormalisation approach to invariant circles in area-preserving maps. Physica D, 7:283-300, 1983. 
http://dx.doi.org/10.1016/0167-2789(83)90131-8 

[Mac89] R.S. MacKay. Flux over a saddle. Phys. Lett. A, 145(8-9):425-427, 1989. http://dx.doi.org/10.1016/ 

0375-9601(90)90306-9 

[Mac94] R.S. MacKay. Transport in 3D volume-preserving flows. J. Nonlinear Science , 4:329-354, 1994. 

[Mat82] J.N. Mather. Existence of quasi-periodic orbits for twist homeomorphisms of the annulus. Topology, 21:457-467, 

1982. http://dx.doi.org/10.1016/0040-9383(82)90023-4 

[MBK + 12] M. Michler, A. Backer, R. Ketzmerick, H.-J. Stockmann, and S. Tomsovic. Universal Quantum localizing tran¬ 
sition of a partial barrier in a chaotic sea. Phys. Rev. Lett., 109(23):234101, 2012. http://dx.doi.org/10.1103/ 
PhysRevLett.109.234101 

[MBR14] C. Manchein, M.W. Beims, and J.M. Rost. Characterizing weak chaos in nonintegrable Hamiltonian systems: The 
fundamental role of stickiness and initial conditions. Physica A, 400:186-193, 2014. http://dx.doi.Org/10.1016/j. 
physa.2014.01.021 

[MD06] K.A. Mitchell and J.B. Delos. A new topological technique for characterizing homoclinic tangles. Physica D, 
221(2):170-187, 2006. http: //dx. doi . org/10.1016/j .physd. 2006.07.027 
[Mei 86 ] J.D. Meiss. Class renormalization: Islands around islands. Phys. Rev. A, 34:2375-2383, 1986. http://link.aps. 
org/doi/10.1103/PhysRevA.34.2375 

[Mei92] J.D. Meiss. Symplectic maps, variational principles, and transport. Rev. Mod. Phys., 64(3):795-848, 1992. http: 
//dx.doi.org/10.1103/RevModPhys.64.795 

[Mei94] J.D. Meiss. Transient measures for the standard map. Physica D, 74(3 & 4):254-267, 1994. http://dx.doi.org/ 
10.1016/0167-2789(94)90197-X 

[Mei97] J.D. Meiss. Average exit time for volume-preserving maps. Chaos, 7:139-147, 1997. http://dx.doi.org/10.1063/ 
1.166245 

[Mei08] J.D. Meiss. Visual explorations of dynamics: The standard mapping. Pramana, Indian Academy of Sciences, 
70:965-988, 2008. http://www.ias.ac.in/pramana/v70/p965/fulltext.pdf 
[MH01] N. Murray and M. Holman. The role of chaotic resonances in the solar system. Nature, 410(6830):773-779, 2001. 
http://dx.doi.org/10.1038/35071000 

[Mit09] K.A. Mitchell. The topology of nested homoclinic and heteroclinic tangles. Physica D, 238(7):737-763, 2009. 
http://dx.doi.org/10.1016/j.physd.2009.01.004 

[MJ13] J.D. Mireles-James. Quadratic volume-preserving maps: (un)stable manifolds, hyperbolic dynamics, and vortex- 
bubble bifurcations. J. Nonl. Sci., 23(4):585-615, 2013. http://dx.doi.org/10.1007/s00332-012-9162-l 
[MM83] R.S. MacKay and J.D. Meiss. Linear stability of periodic orbits in Lagrangian systems. Phys. Lett. A, 98:92-94, 

1983. http://dx.doi.org/10.1016/0375-9601(83)90735-1 

[MM 86 ] R.S. MacKay and J.D. Meiss. Flux and differences in action for continuous time Hamiltonian systems. J. Phys. A, 
19:L225-L229, 1986. http://dx.doi.Org/10.1088/0305-4470/19/5/002 
[MM 88 ] R.S. MacKay and J.D. Meiss. Relationship between quantum and classical thresholds for nrultiphoton ionization of 
excited atoms. Phys. Rev. A, 37(12):4702-4706, 1988. http://pra.aps.org/abstract/PRA/v37/il2/p4702_l 
[MM92] R.S. MacKay and J.D. Meiss. Cantori for symplectic maps near the anti-integrable limit. Nonlinearity, 5:149-160, 
1992. http://dx.doi.org/10.1088/0951-7715/5/1/006 

[MMP84a] R.S. MacKay, J.D. Meiss, and I.C. Percival. Stochasticity and transport in Hamiltonian systems. Phys. Rev. 

Lett., 52(9):697-700, 1984. http://link.aps.org/abstract/PRL/v52/p697 
[MMP84b] R.S. MacKay J.D. Meiss, and I.C. Percival. Transport in Hamiltonian systems. Physica D, 13:55-81, 1984. 
http://dx.doi.org/10.1016/0167-2789(84)90270-7 

[MMP87] R.S. MacKay, J.D. Meiss, and I.C. Percival. Resonances in area-preserving maps. Physica D, 27(l-2):l-20, 1987. 
http://dx.doi.org/10.1016/0167-2789(87)90002-9 

[MMS89] R.S. MacKay, J.D. Meiss, and J. Stark. Converse KAM theory for symplectic twist maps. Nonlinearity, 2:555-570, 
1989. http://iopscience.iop.org/0951-7715/2/4/004 

[MO 86 ] J.D. Meiss and E. Ott. Markov tree model of transport in area-preserving maps. Physica D, 20:387-402, 1986. 
http://dx.doi.org/10.1016/0167-2789(86)90041-2 

[MR14] T. Manos and M. Robnik. Survey on the role of accelerator modes for anomalous diffusion: The case of the standard 
map. Phys. Rev. E, 89(2):022905, 2014. http://link.aps.org/doi/10.1103/PhysRevE.89.022905 
[MS14] R.S. MacKay and D.C. Strub. Bifurcations of transition states: Morse bifurcations. Nonlinearity, 27(5):859, 2014. 
http://dx.doi.org/10.1088/0951-7715/27/5/859 

[MSM13] B.A. Mosovsky, M.F.M. Speetjens, and J.D. Meiss. Finite-time transport in volume-preserving flows. Phys. Rev. 

Lett., 110(21):214101, 2013. http://dx.doi.org/10.1103/PhysRevLett.110.214101 
[Per79] I.C. Percival. Variational principles for invariant tori and cantori. In M. Month and J.C. Herrera, editors, Nonlinear 
Dynamics and the Beam-Beam Interaction, volume 57, pages 302-310. Amer. Inst. Phys., New York, 1979. http: 
//dx.doi.org/10.1063/1.32113 

[P 0 I 88 ] L.V. Polterovich. On transport in dynamical systems. Russ. Math. Sur., 43(l):251-252, 1988. http://dx.doi.org/ 
10.1070/RM1988v043n01ABEH001689 



27 


[PosOl] J. Poschel. A lecture on the classical KAM theorem. In Smooth Ergodic Theory and Its Applications (Seattle, Wa 
1999), volume 69 of Proc. Sympos. Pure Math., pages 707-732. Amer. Math. Soc., Providence, 2001. 

[PP79] P. Pechukas and E. Poliak. Classical transition state theory is exact if the transition state is unique. J. Chem. Phys., 
71:2062-2068, 1979. http://dx.doi.Org/10.1063/l.438575 

[RK94] V. Rom-Kedar. Homoclinic tangles-classification and applications. Nonlinearity, 7(2):441-473, 1994. http://dx. 
doi.org/10.1088/0951-7715/7/2/008 

[RKLW90] V. Rom-Kedar, A. Leonard, and S. Wiggins. An analytical study of transport, mixing, and chaos in an unsteady 
vortical flow. J. Fluid Mechanics, 214:347-394, 1990. http://dx.doi.org/10.1017/S0022112090000167 
[RKW90] V. Rom-Kedar and S. Wiggins. Transport in two-dimensional maps. Arc. Rational Mech. Anal, 109(3):239-298, 

1990. http://dx.doi.org/10.1007/BF00375090 

[RKZ98] V. Rom-Kedar and G.M. Zaslavsky. Islands of accelerator modes and homoclinic tangles. Chaos, 9(3):697-705, 
1998. http://dx.doi.Org/10.1063/l.166444 

[RLBK14] M. Richter, S. Lange, A. Backer, and R. Ketzmerick. Visualization and comparison of classical structures and 
quantum states of four-dimensional maps. Phys. Rev. E, 89(2):022902, 2014. http://dx.doi.org/10.1103/PhysRevE. 
89.022902 

[SLV05] J.D. Szezech, S.R. Lopes, and R.L. Viana. Finite-time Lyapunov spectrum for chaotic orbits of non-integrable 
Hamiltonian systems. Phys. Lett. A, 335(5-6), 2005. http://dx.doi.Org/10.1016/j.physleta.2004.12.058 
[SMOW08] R. Sturman, S.W. Meier, J.M. Ottino, and S. Wiggins. Linked twist map formalism in two and three dimensions 
applied to mixing in tumbled granular flows. J. Fluid Mech., 602:129-174, 2008. http://dx.doi.org/10.1017/ 
S002211200800075X 

[SS87] L.A. Smith and E.A. Spiegel. Strange accumulators. In J.R. Buc.hler and H. Eichhorn, editors, Chaotic Phenomena 
in Astrophysics, volume 497 of Annals of the New York Academy of Sciences, pages 61-65. New York Academy of 
Sciences, New York, 1987. 

[TJDF85] M. Toller, G. Jacucci, G. Delorenzi, and C.P. Flynn. Theory of classical diffusion jumps in solids. Phys. Rev., 
32B (4): 2082-2095, 1985. http: //dx. doi. org/10.1103/PhysRevB. 32.2082 
[Tom96] S. Tompaidis. Approximation of invariant surfaces by periodic orbits in high-dimensional maps, some rigorous 
results. Experiment. Math., 5:197-209, 1996. http://dx.doi.org/10.1080/10586458.1996.10504588 
[UF85] D.K. Umberger and J.D. Farmer. Fat fractals on the energy surface. Phys. Rev. Lett., 55:661-664, 1985. http: 
//dx.doi.org/10.1103/PhysRevLett.55.661 

[Ven08] R. Venegeroles. Calculation of superdiffusion for the Chirikov-Taylor model. Phys. Rev. Lett., 101(5):054102, 2008. 
http://dx.doi.org/10.1103/PhysRevE.77.027201 

[Ven09] R. Venegeroles. Universality of algebraic laws in Hamiltonian systems. Phys. Rev. Lett., 102:064101, 2009. http: 
//dx.doi. org/10.1103/PhysRevLett.102.064101 

[WH91] J. Wisdom and M. Holman. Symplectic maps for the n-body problem. Astronomical Journal, 102(4):1528-1538, 

1991. http://dx.doi.org/10.1086/115978 

[WHK03] M. Weiss, L. Hufnagel, and R. Ketzmerick. Can simple renormalization theories describe the trapping of chaotic 
trajectories in mixed systems? Phys. Rev. E, 67(4):046209, 2003. http://dx.doi.org/10.1103/PhysRevE.67.046209 
[Wig37] E. Wigner. Calculation of the rate of elementary association reactions. J. Chem. Phys., 5:720-725, 1937. http: 
//dx.doi.org/10.1063/1.1750107 

[Wig92] S. Wiggins. Chaotic Transport in Dynamical Systems, volume 2 of Interdisciplinary Applied Mathematics. Springer- 
Verlag, New York, 1992. 

[WilOl] E.J.N. Wilson. An Introduction to Particle Accelerators. Oxford University Press, Oxford, 2001. 

[WM97] J.M. Wendlandt and J.E. Marsden. Mechanical integrators derived from a discrete variational principle. Physica 
D, 106(3-4):223-246, 1997. http://dx. doi. org/10.1016/S0167-2789(97)00051-1 
[WO04] S. Wiggins and J.M. Ottino. Foundations of chaotic mixing. Phil. Trans. R. Soc. Lond., 362(1818):937-970, 2004. 
http://dx.doi.org/10.1098/rsta.2003.1356 

[WW10] H. Waalkens and S. Wiggins. Geometrical models of the phase space structures governing reaction dynamics. Reg. 

and Chaotic Dyn., 15:1-39, 2010. http://dx.doi.org/10.1134/S1560354710010016 
[Zas85] G. M. Zaslavsky. Chaos in Dynamic Systems. Harwood Academic, New York, 1985. 

[Zas02] G.M. Zaslavsky. Chaos, fractional kinetics, and anomalous transport. Physics Reports, 371(6):461-580, 2002. 
http://dx.doi.org/10.1016/S0370-1573(02)00331-9 

[ZEN97] G.M. Zaslavsky, M. Edelman, and B.A. Niyazov. Self-similarity, renormalization, and phase space nonuniformity 
of Hamiltonian chaotic dynamics. Chaos, 7:159-181, 1997. http://dx.doi.Org/10.1063/l.166252 
[ZT91] G.M. Zaslavsky and M.K. Tippett. Connection between recurrence-time statistics and anomalous transport. Phys. 
Rev. Lett., 67(23):3251-3254, 1991. http://dx.doi.org/10.1103/PhysRevLett.67.3251 



